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Abstract Helioseismology has taught us a great deal about the stratification and 
kinematics of the solar interior, sufficient for us to embark upon dynamical studies more 
detailed than have been possible before. The most sophisticated studies to date have 
been the very impressive numerical simulations of the convection zone, from which, 
especially in recent years, a great deal has been learnt. Those simulations, and the 
seismological evidence with which they are being confronted, are reviewed elsewhere 
in this volume. Our understanding of the global dynamics of the radiative interior of 
the Sun is in a much more primitive state. Nevertheless, some progress has been made, 
and seismological inference has provided us with evidence of more to come. Some of 
that I summarize here, mentioning in passing hints that are pointing the way to the 
future. 
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1 Introduction 

The most studied dynamics of the Sun occurs in the convective envelope, where the 
timescales are human. The large-scale convective flow and its associated angular- 
momentum transport is to some degree accessible to seismological probing. However, I 
sh all pay little atten tion to that here, because i t is reviewed elsewhere in this volume 
bv iBrun et al.l (l2015h and lHanasoee et al.l (l2015h . who discuss principally the results of 
the impressive modern simulations and the physics that has been gleaned from them. 
Instead, I shall try to reflect on the global behaviour of the Sun, recalling some of the 
issues that have concerned heliophysicists in the past - issues that to some degree have 
been thought by commentators not in the field to have been resolved, yet in the light of 
our increasingly sophisticated helioseismological findings have left some room for more 
than a modicum of uncertainty. 

In order to establish the what-one-might-call standard view of the Sun, I shall first 
describe briefly what we believe to have been the principal dynamical processes that 
have influenced the evolution of the Sun to its present state. Only then do we have 
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a basis for discussing the extent to which those beliefs should really be relied upon. 
Quite naturally, I first adopt the simplest of descriptions, consistent with our general 
understanding of stellar evolution; that seems to provide a gratifyingly accurate first- 
order class of models. There have been attempts to standardize some versions of those 
models, principally by John Bahcall. However, a universal standard has never been 
produced, because heliophysicists have not agreed on the manner in which the physics 
should be sim plified. The discord has been a ired many times, such as during lAU Col- 
loquium 121 (IBerthomieu and Cribierl Il990l i . at which Evry Schatzman (reported by 
lOoughl [l990bh favoured the inclusion, by some means or other, of all the pertinent 
generally accepted physics, including macroscopic dynamics (represented in some pre¬ 
cisely stated, albeit approximate, manner), in order to obtain the most realistic model 
possible. By contrast, John argued for extreme dynamical simplicity, ignoring all the 
effects of macroscopic motion except the heat flux by thermal convection (yet exclud¬ 
ing the momentum and kinetic-energy fluxes), parametrized by a local mixing-length 
formula, in order to produce a very straightforward model that can be used as a stan¬ 
dard benchmark, although in practi ce that standard has evolved with t ime, even for 
John. Here, I shall adopt Model S of IChristensen-Dalsgaard et al.l JlflQfil ') as my refer¬ 
ence, because the assumptions upon which it was built have been stated clearly, and 
because, unlike many other solar models, it has been produced with sufficient care for 
reasonably precise, well defined (adiabatic) oscillations to be computed from it. 

I do not go into the details of the early, pre-main-sequence, dynamics. That remains 
a complicated area of (active) research which continues to benefit from advances in 
(principally numerical) simulation. Suffice it to say that, because the age of the Sun is 
very much greater than any of the dynamical timescales, its present structure is hardly 
likely to have any recall of the precise manner in which the Sun condensed from the 
interstellar medium, save for the relic chemical composition, and the total mass and 
angular momentum. It is normally presumed that the chemical composition was either 
initially uniform, or that mixing during Hayashi contraction homogenized any initial 
internal inhomogeneity well before the star settled onto the main sequence. However, 
that assumption is open to question. It is also normally presumed that there was no 
substantial mass loss nor accretion, although that presumption is also questioned from 
time to time, with good reason. It is generally believed that, in common with other 
stars, there has been some loss of angular momentum, particularly in the immediately 
pre-main-sequence and early main-sequence stages of evolution (and even more before), 
as is evinced by the decreas e with age of the rotational velocities of the photospheres of 
stars in young clusters ('e.g. lSkumanic 3ll9z3). However, that process is unlikely to have 
had a great effect on the state of the Sun today (in other words, it is of little concern 
that the angular momentum in the very early days might have been rather greater than 
it is today), and is accordingly rarely discussed. Nevertheless, a word or two about it 
later in my discussion may not be out of place. It may have been responsible for a 
degree of mixing of the products of nuclear reactions in the core, especially in early 
times. Be that as it may, the essentially static effect of rotation, via the centrifugal 
force, is very much smaller than the pressure gradient and the force of gravity (as 
also are Maxwell stresses and the momentum flux due to large-scale convection, except 
in the very outer, near-surface, layers) so to a first approximation the Sun can well 
be regarded as being spherically symmetric. Oblateness of the figure resulting from 
rotation can be treated subsequently as a small perturbation, as I discuss later. So can 
the effect of mass loss. 
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2 Standard main-sequence evolution 


The basic principle behind the evolution is quite straightforward and very robust. 
Nuclear reactions in the core convert hydrogen into helium, reducing the number of 
particles - nuclei, ions and electrons - per unit mass of the stellar material, reducing 
the pressure at a given temperature, permitting the star to contract under gravity 
and thereby raising the temperature in the innermost regions to maintain hydrostatic 
balance. The temperature- and density-sensitive nuclear reaction rates are augmented 
with this contraction by more than the reduction due to the loss of hydrogen fuel, 
essential to maintain global stability, so the luminosity of the star increases. This 
basic result is an inevitable consequence of thermonuclear physics and the conditions 
for hydrostatic balance, and is beyond challenge (provided, of course, that standard 
nuclear and gravitational physics are accepted). 

The theoretical variation of the luminosity L with time t is given approximately by 


L(t)~ 






( 1 ) 


( Goughl[l977bh with p = L{0 )/Lq — 1 ~ 2/5, where the origin of time t is at the start 
of the evolution on the main sequ ence. It is robust against minor variations in the 
assumptions behind the model (e.g.|Gough||l98S ;lBahcall et al.l[20oih. A ‘derivation’ of 


this formula is presented bv iGoughI 1 19881 . 1990c ) based on homology scaling, although 


the constant 2/5 was obtained by adjusting an initial theoretical estimate to account for 
the deviation from homologous variation due to nuclear transmutation such as to render 
equation (1) consistent with numerical stellar-evolution computation, so the relation 
should perhaps be regarded more as a physically motivated interpolation formula. The 
formula has been generalized to take account of a putative temporal variation of the 
gravitatio nal constant, G , and has also been adapted to accommodate main-sequence 
mass loss I Goughl[l990cl l. thereby summarizing the numerical computations that have 
been carried out both before and since. 

One of the reasons for considering a temporal variation in the gravitational con¬ 
stant was to obviate what has bee n called the faint-Sun pr oblem (some have even 
called it a paradox). It was posed bv ISagan and MullerJ ( 1972l ~). who argued from sim¬ 
ple equilibrium energy balance that early in the Sun’s main-sequence evolution the 
Earth would have been completely glaciated had the Sun evolved essentially accord¬ 
ing to equation ©, which is cou nter to geolo gical finding. The solar lum inosity is a 
steeply increasing function of G I Teller! Il948l '): L oc ® (lGamowlll967t l. Therefore 
adopting an appropriate small decline of G with time can almost annul the rise in L 
described by equation ([T]): if, for example, G{t')/G{tu) = {t'/tu)~'^ (where t' is time 
measured from the Big Bang and tu is the age of the Universe), then the irradiance 
on Earth would have been very nea rly constanln if —G/ G\t'=t.. — 1.2 x 10^^ y, almost 
irrespective of the value of q (e.g. iGoughI[19881 . Il990dl . However, this is not strong 
evidence in favour of there being a temporal variation in G, because it is unlikely 
that the climatological equilibrium energy-balance assumption is correct. Apparently 
more sophisticated meteorological arguments have even been perceived to exacerbate 
the issue. For example, in the summer of 1973, motivated by a prediction that the 
Sun’s lumi nosity might have decli ned suddenly by a few per cent some million or so 
years ago ( Dilke and Goughlll972[ l. Tzvi Gal-Ghen and I (unpublished) engineered a 
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having deviated from its initial value by no more than 0.5 per cent. 
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5-per-cent reduction in solar irradiance in the Global C irculation Model (GCM) of 
the Earth’s atmosphere i Kasahara and Washingtonlll97lh at the US National Center 
for Atmospheric Research. As expected, the Earth became completely glaciated; but 
what was not expected is that when the irradiance was subsequently restored to its 
present value, the ice that had just been produced would not melt. That result is ev¬ 
idently contrary to observation, because the Earth is not completely glaciated today. 
The moral, for me at least, was to confirm that one must always beware of the products 
of complicated computer programmes, particularly when they are state of the art yet 
depend, as does the GCM, on rather primitive physics. The many physical processes 
that the GCM had to account for were ill understood - as indeed many are still - and 
their influence on the atmospheric dynamics had to be parametrized; the controlling 
parameters whose values were not previously known were measured, where possible, 
and then held constant in the model. That was the most obvious flaw. We know that 
the stability and evolution of any dynamical system can depend quite sensitively, and 
sometimes critically, on the nature of the constraints that are imposed; evidently a 
meteorological model, however sophisticated, with processes whose controlling param¬ 
eters do not change appreciably on a timescale of days or weeks could not possibly be 
expected to apply to the climate on timescales of millions of years. Indeed, Dilke and 
Gough had presumed that the response of the Earth’s climate to irradiance variation 
was quite different on a timescale of 10®y from that on a timescale of 10®y. The pro¬ 
cesses responsible for s uch timescale dependen c e hav e been studied subsequently by , 
amo ngst many others, M argulis and LovelockI ( 1974h . iLovelock and Margulij (Il974l ') 
and ILovelock and Watson 1 1982h . who likened the long-term stability of the climate 
to the work of Gaia, the ancient Greek goddess of Earth. The cause of the suspected 
sudden relatively recent decline in L was itself a dynamical process, considered to have 
resulted from a nonlinear instability of g modes in the Sun’s core, the theory of whose 
development also depending critically on what is held constant, a matter to which I 
shall return later. 

The computation of the structure and evolution of (standard) solar models is 
straightforward once the equation of state, the nuclear reaction rates and the opac¬ 
ity of the stellar material have been specified. There is also a matter of determining 
the mean stratihcation of the convection zone: to that end some form of mixing-length 
formalism is normally used to specify the energy (usually just heat) flux and, some¬ 
times, the Reynolds stres^. For most purposes the deficiencies in such a procedure are 
of little concern to the overall structure of the star, because throughout almost all of 
the convection zone the stratification is adiabatic and the Reynolds stress is negligible; 
it is only in the outer boundary layer where the details of the formalism matter. In 
practice the mixing-length formalism contains at least one explicitly adjustable con¬ 
stant which can be calibrated to determine the value of the adiabatic constant deep 
in the convection zone that is required to reproduce the observed radius of the Sun; 
calibrated in that way, the outcome is a function of chemical composition, and it for¬ 
mally determines the depth of the convection zone. Of course, that is not to say that 
in reality the dynamics of the upper convective boundary layer depends explicitly on 
the chemical composition in the manner determined by the solar calibration. 

The chemical composition is normally specihed by two of three parameters: the 
relative abundances X, Y, Z by mass of H, "^He and of all other elements combinec|3; 


^ yet hardly ever the flux of kinetic energy 

® usually referred to as ‘heavy elements’, even though they include ®He 
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Fig. 1 Square of the sound speed, c^, and temperature, T, in two solar models computed 
with a local mixing-length formalism and having differe nt initial heavy-element abundance s 
Zq: one (continuous curves) is the standard Models S of IChristensen-Dalsgaard et al.l (Il996lh 
having Zq = 0.020 and a corresponding initial helium abundance Yq = 0.27; the other (short- 
dashed curves) has Zq = 0.001, Yq = 0.16 and has been continuously contaminated with 
heavy elements at the the surface at such a constant rate as to have a surface abundance 
Zs = 0.020 today dChristensen-Dalsgaard et al.lll979alll979bll . The depression in (relative 
to T) in the central regions of the cores of the evolved models is a result of the augmentation 
of the mean molecular mass in the core produced by the transmutation of hydrogen into 
helium. Included, for comparison, as a long-dashed curve, is the square of the sound speed of 
the higher-Y model at zero main-sequence age, which, of course, has no such depression. The 
vertical arrows (with line styles corresponding to the models to which they refer) mark the 
bases of the convection zones of the two present-day models, where the second derivatives of 
and T are discontinuous. The square of the sound speed in the Sun, inferred seismologically, 
is drawn also as a continuous curve, not quite reaching the centre r/R = 0 of the star; it is 
barely distinguishable from in the higher-Y model (after [G^ iacni)). 
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they satisfy X-\-Y -\-Z — 1. (The relative abundances of the elements incorporated into 
Z are provi ded principally, but not completely, by spectroscopic analyses of the Sun’s 
atmosphere l Asplund et aL|[2009l : ICaffau et al.ll2009h .') Stellar evolution theory provides 
a further relation between them from the requirement that the current luminosity Lito) 
agrees with the value measured (assuming the age Iq to be known). That leaves a single 
infinity of apparently acceptable models, which can be labelled by any one of X, Y 
and Z. I record that along the sequence of these models, both Y and the depth dc of 
the convection zone increase monotonically with Z. 

Before helioseismology, there was just one measured quantity associated with the 
Sun that in principle could be used to select the most appropriate model, namely the 
neutrino luminosity Lu (usually expressed as an associated flux Fu at lAU computed 
under the assumption that neutrinos do not decay, which itself is expressed as a cap¬ 
ture rate by whatever neutrino detector is being considered). Theoretically, Lv is a 
monotonic increasing function of Z. Therefore a unique model can thereby be chosen. 
However, that model left many people uncomfortable, and was deemed unacceptable by 
most astronomers and astrophysicists because the value of Y in that model was much 
lower than that observed in open stellar clusters containing apparently Sun-like stars; 
moreover, it was lower than the amount of helium created in generally accepted mod¬ 
els of Big-Bang nucleosynthesis. Astronomical opinion prevailed, and the discrepancy 
posed what was accordingly named the solar neutrino problerrQ. 

Perhaps the first genuinely seismological inference to be drawn about the st ructure 
of the Sun was from a crude calibration of the depth dc of the convection zone I GoughI 
Il977al l. It implied a value of Y even greater than those preferred at the time, exac¬ 
erbating the neutrino pro blem. The infe r ence w as soon confi rmed with more detailed , 
nu merical, calculations by [Rhodes et al.l ( 1977h . and later bv iBerthomieu et al.l ( 198fll ') 
an d lLubow et al.l (Il980lj . who addressed the sensitivity of the eigenfrequencies to uncer¬ 
tainties in the background model envelope. Subsequently, a helioseismological inference 
of th e sound speed throughout most of the solar interior I Christensen-Dalsgaard et al.l 
119851 1 not only provided a consistent direct measurement of dc but also ruled out the 
Tiy-reproducing low-T model. A more recent, and cleaner, demonstration of that in¬ 
ference is reproduced in Figure 1. It convincingly demonstrated that the resolution of 
the solar neutrino problem must lie in nuclear or particle physics, and not in stellar 
physics, although that conclusion was not accepted immediately by most of the he- 
liophysics community w ho did not yet appreciate the power of helioseismology (e.g. 
iBahcall and Ulrichlll988h . I shall say no more ab out that here because i t is on ly indi¬ 
rectly associated with interior dynamics (but see lCummine and Haxtt^ 1 199(]l l.l 


3 Angular velocity and the solar oblateness 


The most obvious indicator of global solar dynamics is the angular velocity f2(r, 6, t) 
(I adopt spherical pol ar coordinate s (r,^ ,(j>)). A seismological determination of Q in 
the Sun’s envelop e by Schou et al. l lll998ll is illustrated in Figure 4, and a schematic 
representation by IChanlin et al.l (Il999l l extending to the centre is plotted in Figure 
2. Roughly speaking, the angular velocity is independent of radius and increases with 
colatitude 9 in the convective envelope (r > re), and it is approximately uniform in 
the radiative envelope beneath, except possibly in the energy-generating core (r < rc) 


^ rather than, for example, the solar abundance problem 
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where there is (slig ht) evidencj^ that the spherically averaged 17 is lower than 17 in the 
radiative envelope (lElsworth et alJIlQQRh . The transition at th e base of the co nvection 
zone, called the tachocline (Spiegel and Zahn, 1992; see also ISniegell 1 1972h ). is too 
abrupt to be properly resolved seismologically, as is the transition, if there is one, at 
the edge of the core. 

It is common practice to expand the latitudinal dependence of 17 in even powers of 
jj, — cos9-. 

^ 2k 


n(r,e,t) = ^ 17 fe(r,t)/i 


( 2 ) 


fc >0 


where, on the whole, varies weakly with r and t, except in the tachoclin^E- It 
is interesting and, so far as I am aware, unexplained, that for most values of r the 
magnitude of 17^ is particularly small for all k > 2. The coefficients of the terms of 
lower degree are observed to vary somewhat with the sunspot cycle, but that variation 
has been convincingly detected only in the convection zone, so I refrain from discussing 
it in a ny detail here. There is also a variation with a characteristic period of about 2 
years ( Broomhall et al which has been potentially misleadingly called the quasi¬ 

biennial oscillation - it can hardly be compared with the terrestrial disturbance with 
the same name, having the acronym QBO, because the latter is driven by differential 
gravity-wave dissipation, a process which I explain in §10, and which can hardly be 
sustained in the convection zone of the Sun, where gravity waves cannot propagate. 

The first conclusion to be drawn from the earliest determination of 17 was an 


estim ate of the quadrupole moment J 2 of the external gravitational field (iDuvall et al.l 
Il984l ). It is induced by the centrifugal force acting on the stellar material, causing the 
mass distribution to be oblate. Therefore J2 can be represented as an appropr iately 
weighted average of 17^ over the volume of the star ( Goughlll98ll : |Piipersl|l998h . the 
weight function F being determined by linear perturbation about the nonrotating state. 
The radial dependence of the component of F pertaining to the spherically averaged 
17 is plotted in Figure 3. Knowledge of J2 is essential for testing theories of gravity 
from measurements of planetary precession and spacecraft orbits, for non-Newtonian 
theories produce a perturbation to Newtonian orbits that at present are observationally 
indistinguishable from that which is produced by a deviation from spherical symmetry 

induced by an appropriate internal centrifugal fo rce. _ _ 

Ea rly a ttemp t s to measure J2 were made by lAuwerj ( 189ll ) , lAmbronn and Schurl 
( I9O5I ) and IPoorl 1 1905alll905lJ) from direct observations of the presumedly centrifu- 
gally induced visual oblateness, defined as Av = {Re — Rp)/R, where Re, Rp and R 
are respectively the equatorial and polar radii and their avera gjlt in retrospect those 
observations were inconclusive. More modern measurements bv lpicke and Goldenber3 
1 19671 . I1974I ) suggested that Av is about 5 times greater than what one would ex- 
pect from a distortion by a uniform rotation of the Sun’s interior, consistent with the 
iBrans and Dicke dlMlI) theory of grav ity, appropriately calibrated. However, that was 
challenged by Hill and Stebbin'i ( IQTSl ) , whose measured value was consistent with uni¬ 
form rotation. These and subsequent measurements have been reviewed by, for example, 


^ which is commonly doubted 

® Often, for computationa l convenience, f? is expanded in orthogonal polynomials, such as 
Clebsch-Gordon coefficients dRitzwoller and Lavelvlfl99ffi 
^ The oblateness Ay depends also on J 4 and the higher moments, but the additional contri¬ 
butions appear to be less than the observational uncertainty, so for clarity I do not take them 
explicitly into account here. 
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Fig. 2 Angular velocity through the Sun, plotted against radius at three different latitudes. 
The continuous curves are the expectations (averages), and are flanked by dashed curves with 
added or subtracted standard errors propagated from the estimated errors in the data. The 
tachocline is not well resolved; nevertheless its base appears to be spherical, as one would 
expect, and penetration of the shear into the convection zone preferentially at high latitudes, 
rendering the overall shear layer prolate, is discernible. The angular velocity in the radiative 
envelope (outside the energy-generating core) is essentially unifor m, within the uncer tainty, 
and there is a hint that the core is rotating more slowly ffrom IChaplin et ahlflOQ^ . The 
confluence at r/i? ~ 0.25 of the values of Q at the three different latitudes and the almost linear 
dependence beneath are both consequences of the regularization, which, because the seismic 
data are too scant to resolve the angular velocity well, dominates the inversion procedure. All 
that can be inferred is that an equatorially biassed latitudinal average of the angular velocity 
in the core appears l ikely to be lower than the angular velocity of the radiative envelope (cf. 
lElsworth et al.lll995l^ : therefore the apparent uniformity of Q with respect to latitude in the 
core cannot be regarded as evidence for an absence of latitudinal variation. 


iRozelot et al.l (|2Qllh . The shape measurements were very difficult, partly because Z\v 
is very small, of the order of the centrifugal parameter A = /GM ~ 2 x 10“^, 

where M is the mass of the Suu, and partly because we now know that the oblateness 
ziv of the visible solar disc is dominated not by the oblateness of the gravitational 
field caused by the action of the centrifugal force on the dense interior material, but 
principally by the oblateness Aq caused by the direct effect of the centrifugal force 
on the diffuse visible surface layers. Therefore A,p would need to be determined as the 
relatively small difference between A-v and Aq , which is an intrinsically unreliable pro¬ 
cedure. The difficulties in the entire measurement procedure have been compounded 
by the fact that the relation between the oblateness of the observed radiative inten¬ 
sity and the oblateness of, say, the surfaces of constant pressure in the photosphere is 
contaminated with brightness variations due to sunspots and magnetically generated 
excess emission, such as faculae. Accounting for t hese is an uncertain process, as is 
evident from even the most recent investigations i Lefebvre et al ] I2OO7I : iFivian et al.l 
l2008l:lKuhn et al.ll2012h . 




















9 



Fig. 3 Superposed on the setting Sun is a plot of the weighting function F(x) in the approx¬ 
imate formula J 2 — J -Fi? da;, where x = r/R and i? is the spherical average of the (lati- 
tudina lly mildly varying) internal angular velocity of the Sun, in units of s“^ (after [G^ iE 
l2012aD . The function F is small near the centre of the Sun, where the ratio of centrifugal force 
to angular velocity is low, and near the surface where the low density, however oblate, makes 
only a minor imprint on the gravitational potential. 


The advent of helioseismology changed the situation dramatically, because the seis¬ 
mic signature pertinent to determining J2 can be measured much more accurately: that 
signature is the rotational frequency splitting, whose magnitude is of order 2mf2, where 
m is the magnitude of the azimuthal order of a seismic mode; modern measuremen t 
error of mean dipole-mode splitting, for example, is approximately 1% 1 Howel [2OO9I I . 
which offers some idea of the precision to which 17 can be inferred. Furthermore, the 
dynamical relation between 17, which is inferred directly from the frequency splitting, 
and J 2 is not influenced significantly by non-dynamical variables such as radiative 
i ntensity. From eve n the first primitive measurement of the interior angular velocity 
( Duvall et al.l[l983l it was evident that the prediction of the non-Newtonian component 
of the precession of the perihelion of the orbit of Mercury is con sistent with Genera l 
Relativity. Subsequent measurements, such as those reported bv ISchou et al ] dlQQSll . 
tightened that conclusion, and will in future provide more stringent tests of theories of 
gravitation once finer orbital measurements become available. 


It should be noted that the only unambiguous seismic signature of rotation comes 
from inertial (such as Coriolis) effects on the nonaxisymmetric (m > 0) modes, which 
are tiny in the evanescent region of those modes: namely, near the axis and, doubly, 
near the centre (because the degree I cannot be less than m), which explains the white 
region in Figure 4. That region coincides, accidently, with one of the regions in which 
the oblateness kernel (whose spherical-average component is illustrated in Figure 3) is 
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Fig. 4 Colour rendering of the rotation rate of a meridional quadrant of the Sun. The scale 
is in nHz. The continuous quarter circle is located in the photosphere, the dashed quarter 
circle at the base of the convection zone. The white area denotes the region in which the 
rotational splitting data provide only scant indication of the angular velocity (from the MDI 
I mage Gallery: http:/ /soi.stanford.edu/results/Triana/modslice.ps constructed from the work 
of lSchou et al.l ||199^ 1 V 


small, because the centrifugal force is small. Therefore the principal uncertainty in the 
inference of O contributes little to the uncertainty in the value of that is derived. 


4 A dynamical issue raised by recent oblateness measurements 

Taken at face value the most recent direct measurements raise an interesting dynamical 
question. After doing their best to remove facular contamination by rejecting exces¬ 
sively bright regions at the solar limb, which are equatorially concent rated and therefore 
contribute positively to the oblateness of the intensity distribution, iKuhn et al.l (l2012l'> 
found a residual visual oblateness that is actually less than Aq by about 8 x 10”^. The 
observations by Fivian et al. (2008) led to a similar, although less extreme, resul10. 
So is the Sun gravitationally prolate? The only way that that could be is probably 
for the radiative interior to be constricted by a predominantly toroidal belt of mag¬ 
netic field (with an associated poloidal component, for stability) encircling the rotation 

® In their paper, iFivian et al.l ll2008f) actually claim a value that exceeds Aq by about 
2 X 10“^, which they offer as an estimate of but their result was obtained with Dicke’s 
(1970) outdated underestimate of Aq. 
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axis, of intensity that would need to be lO^T or more. Is that possible? What are the 
alternatives? 

Of course there is always the possibility that the observations were inadequately 
analysed. Kuhn et al. explicitly rejected bright facular regions. Did they take adequate 
care to reject darkened, presumably penumbral regions whose inappropriate inclusion 
could have rendered the radiative intensity more prolate? Fivian et al. accounted for 
faculae by using a proxy indicator for rejecting data and extrapolating to 100% 
rejection; there is a fear in some quarters that the resulting statistics are inadequate for 
rendering the outcome reliable. A naive theorist insufficiently familiar with the details 
of the observations might simply take the difference between the two inferred values of 
the uncontaminated visual oblateness, about 8 x 10”^, as a plausible estimate of the 
uncertainty, and, noticing that it is nearly three times the value of A,/, obtained from 
helioseismology, set the matter aside. 

It is nonetheless interesting to enterta i n the idea that there is a prolate contribution 
to the distortion of the surface. As IPic y lll97nll has emphasized, that must necessarily 
occur in the seen photospheric layers. A solar physicist’s immediate reaction is likely to 
be to invoke a superficial magnetic field. The photospheric mean (spatially smoothed) 
field appears to be principally poloidal and dipolar at most times, and of insufficient 
intensity to maintain the pole-equator photospheric radi us difference of the order of the 
that would be required to explain the findings of iKuhn et al.l ( 2012l ~l. Moreover, 
it would enhance, not reduce, the visual oblateness. On the other hand, a preferential 
suppression of polar convection by only a mere 10”^ per cent, which is arguably more 
plausible, would be sufficient to induce the required reduction in oblateness if the local 
static adjustment of the star in only the vertical direction were to matter. What would 
be the dynamical consequence? The poles would be elevated relative to a surface that 
is everywhere perpendicular to the combined gravitational and centrifugal forces, and 
if convective Reynolds stresses acted like a scalar viscosity, which is not uncommonly 
assumed yet which is almost certainly incorrect, photospheric matter would flow down¬ 
hill towards the equator, which is contrary to the direction observed. The dynamical 
problem is evidently not straightforward. Yet it is an interesting and evidently im¬ 
portant problem, for, irrespective of the oblateness issue, a cogent explanation of the 
poleward meridional flow in t he outermost layers of the convection zone is lacking (e.g. 
iToomre and ThomosonlboiiTh . 


5 Spin-down 


The solar wind is rotationally coupled to the Sun via a large-scale magnetic field. It is 
caused to rotate roughly at the solar rate out to about 5 solar radii, thereby removing 
angular momentum from the Sun; beyond, angular momentum in the wind is more-or- 
less conserved. Thus the rotation of the outer layers of the Sun must be decelerating, in 
common with inferences for other stars drawn from observations of rot ational spectral- 
line broadening in open clusters of difference ages ( Skumanichlll973l . An interesting 
dynamical question that arises from this process is how far into the Sun this deceleration 
penetrates. Of course we now know the answer from seismological measurement, but 

the dynamical issues, to some extent, remain. _ 

The matter was debated in the late 1960s after ( Dickel[l964h had tried to maintain 
that the Sun’s core is rotating so rapidly (with a period of a little over a day) as to in¬ 
duce a gravitational oblateness of magnitude sufficient to sever the agreement between 
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the observed residual precession of the orbit of Mercury and the prediction of General 
Relativity. Dicke pointed out that the global viscous diffusion time Tv = F? jv (where v 
is a mean kinematic viscosity) exceeds even the age of the Universe, and that therefore 
the Sun’s core is rotating essentially at its pre-main-sequence rate. He argued that the 
shear could be stable to the Richardson criterion, so shear turbulence wo uld not add to 
the vi scous stress, and that Maxwe l l stres ses could be insignihcant too. iHoward et al.l 
( 1967h and ISretherton and Soiegell 1 1968h pointed out a fundamental flaw in Dicke’s 
claim, using simple dynamical models to demonstrate that meridional advection, and 
not viscous stress, is likely to be the dominant agent transpor ting angular momentum 
through the body of the Sun, in a process called spin-down I Greenspan and HowardI 
Il963l'l . as in a stirred cup of tea. 

Interesting ly, the spin-down process had been discussed qualitatively long before by 
lEinsteinl (192a) as an explanation of the meanders of rivers (unwittingly defendin g his 
theory of General Relativity), and quantitatively bv iBondi and LvttletonI 1 1948h in a 
discussion of the deceleration of the Earth’s core. I describe the broad principles briefly 
here, because, as I shall mention later, the general conclusion may find astrophysical 
application elsewhere. 

Consider a stationary cylindrical vessel of radius R containing water initially ro¬ 
tating (approximately) uniformly w ith angular veloc ity 17 about its vertical axis, a 
situation discussed, for example, bv ICreenspanl ( 1969l'> . The situation is illustrated in 
Figure 5. The fatter (red), horizontal, double arrows represent the centrifugal force, 
which is balanced by a pressure gradient directed away from the axis. The excess pres¬ 
sure Ap near the outside wall of the container supports a greater head of water, and 
therefore the upper, free, surface of the water is concave upwards. Surfaces of constant 
pressure are not normal to gravity, a situation known as baroclinicity. Viscous stresses 
slow the rotation of the water near where it touches the container, particularly at the 
bottom, thereby reducing the centrifugal force. That leaves an unbalanced component 
of the pressure gradient which drives a spiralling inward flow with characteristic speed 
u, indicated in Figure 5 by the thinner (green) arrows, in a thin boundary lay er (of 
characteristic thickness S R), now called an Ekman layer (after lEkma iil ligosli , near 
the bottom of the container. It is there that angular momentum is removed from the 
watei[^. In the boundary layer the Goriolis force (in a frame rotating with the bulk of 
the fluid) is balanced by the viscous stress in the shear: Qu ~ v5~^u. The radial pres¬ 
sure gradient, of order ApjR, is also balanced by viscous stress, which determines the 
boundary-layer velocity: Ap/R ~ RQ^ ~ vujS^. The return flow - upwards and out¬ 
wards in the bulk of the fluid where viscous stresses are negligible, with characteristic 
velocity w ~ u5/R - slows the rotation in the bulk of the fluid mainly by angular- 
momentum conservation, on the timescale r ~ R/w. This spin-down (equilibration) 
time is thus r ~ R?/5u ~ R /\/it is the geometric mean of the global viscous 
diffusion time Tv = R?/v and the characteristic dynamical time 17“^, which is much 
shorter than the viscous time. More accurate formulae for 5 and r can be obtained by 


® Viscous stress operates also on the side walls, but there the boundary layer is not as thin 
as that at the bottom of the container, and removes negligible angular momentum (just as 
Bondi and Lyttleton had found, in the case of spin-down in a sphere, that negligible angular 
momentum is removed near the equator). The bottom boundary layer is thinner as a result of 
the vertical ‘rigidity’ imparted on the fluid immediately above by the vortex stretching (which 
is intimately related to a tendency towards local angular-momentum conservation) produced 
by the shear, and which is also responsible for the better-known Taylor-Proudman theorem 
for steady incompressible inviscid flow. 
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Fig. 5 Spin-down in a stirred cup of tea. The free surface of the fluid is denoted by the curved 
black line, the cylindrical vessel containing the fluid by the vertical and horizontal lines. The 
thicker arrows (red) with black heads indicate the centrifugal force acting on the body of the 
rotating fluid; that force is essentially absent from the Ekman boundary layer near the base 
of the container, which hardly rotates. The thinner (green) arrows indicate the direction of 
the induced meridional flow, which, in the main body of the fluid, is essentially inviscid and 
therefore angular-momentum conserving. Angular momentum is lost from the fluid principally 
by viscous transfer in the thin Ekman boundary l ayer. Such a boundary layer is not essential 
to the global mechanism of spin-down, however: as iBretherton and Soieg^ lll968l') have argued, 
penetration of the meridional flow into the Sun’s convection zone could provide an (even more) 
effective mechanism for divesting angular momentum. That process can be modelled on Earth 
by spreading a layer of heavy beads over the bottom of the teacup. 


a local analysis of the boundary layer, which relates the vertical shear in the angular 
velocity to the upflow velocity out of the Ekman layer I Greenspanlll96^) consequent 
on what is now known as gyroscopic pumping. 

An important observation of the process is that the flow near the bottom of the 
container is inwards, towards the central axis. That is why loose tealeaves in a stirred 
cup of tea migrate to the middle of the cup, and silt at the bottom of a slowly flowing, 
meandering, river is t ransported fro m the outer to the inner bank of a bend, accen¬ 
tuating the meander ( Einsteinlll926l) . It is evident from the balance of forces that if 
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instead the container were caused to rotate more rapidly than the fluid, a flow with 
the same geometry would be induced, but with its direction reversed. 

The principal importance of this process for spin-down is not in the details of the 
Ekman layer, but simply in the fact that the layer extracts angular momentum from 
the fluid, leaving an unbalanced pressure gradient within the layer to drive a large- 
scale flow towards (or away from, were the teacup to be rotating faster than the tea) 
the axis. I have included the account of the boundary layer here, however, simply to 
complete my discussion of some interesting physics. 


Returning now to the discussion of the Sun, one must first beware, quite generally, 
of attributing the properties of oversimplified models to complicated situations with¬ 
out careful consideration of the implications of those simplifications, especially when 
dynamics is involved. Good models might exhibit some of the physical processes in 
operation, but it must always be appreciated that they may be no more than merely 
illustrative. It is therefore not uncommon for arguments based on the understanding 
gleaned from them to be challenged, in disbelief of the extension of the domain of 
applicability necessarily required for addressing the matter in hand. 


For example, IPic ^ lIlQbit) pointed out that the Sun is no cup of tea, and that the 
base of the convection zone does not exert a stress in t he manner of the base ( o r top ) 
of a rigid container to produce a diffusive Ekman layer. iBretherton and Sniegell ( igbSl i 
countered by explaining that an Ekman layer is not essential: the pertinent agent is 
the induced meridional flow which advects angular momentum essentially inviscidly 
throughout the body of the fluid; that flow can penetrate into the convection zone 
where it can divest its angular momentum even more efficiently than in a viscous 
boundary layer. They illustrated the process by modelling the convection zone as a 
rigid porous medium!^ They modelled the radiative zone as an incompressible fluid, 
since the details of how exactly angular momentum is transported in the essentially 
inviscid radiative interior, be it compressible or not, is of relatively minor importance. 


Dicke also tried arguing that the stable stratification of the radiative interior per¬ 
forms the c rucial role of inhibiti ng large-scale, angular-momentum-advecting fl ow - as 
indeed had Howard et al.l (Il967l 'l recognized already - confining it to a shallow iHoltonI 
( 1965alll965bl l layer immediately beneath the convection zone, thereby insulating the 
core from the deceleration of the surface (a conclusion which is v alid only on timescales 
shorte r than the shortest diffusion - here thermal - timescale); [McDonald and Dickel 
(ll9B7l ~) illustrated the process experi mentally, concluding that the stratification of the 
Sun s hould preclude core spin-down. IClark et al.l ( 19691 1 and iModisette and Novotnvl 
(ll969|) concur r ed. A pparently oblivious of the arguments that had been presented by 
[Howard et al.l ([l9B7h . Dicke had failed to recognize that actually the impediment is 
thermally moderated - negative buoyancy generated by vertical adiabatic motion can 
be annulled by thermal diffusion - so one must estimate the buoyan c y ann ihilation time 
in order to assess the validity of the simplified model. [R.oxbur^ il9B4h had already 
pointed out the importance of thermal diffusion, positing that the diffusively moderated 
spin-down time is the (thermally driven) Eddington-Sweet timescale, which is greater 
than the age of the Sun. But Howard, Moore and Spiegel considered, more realistically, 
the thermal control of the (dynamically driven) flow, estimating the spin-down time to 
be at most of order only lO^yr. To fluid dynamicists at the time of the 1960s debate, the 


and subsequently demonstrated it in the laboratory (unpublished) in a rotating beaker of 
water containing several layers of glass beads. 
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balance of the evidence seemed to favour substa ntial global spindownr I but the case 
had not yet been proven. Subsequent analysis bv ISniegel and ZahnI (11992 1 of a similar, 
though not identical, situation relating to the seismologically observed near-uniform 
rotation of the radiative interior in the face of a differentially rotating convection zone, 
to which I turn my attention in the next section, has added substantial dynamical 
support to the evidence. And, of course, the seismological findings themselves negate 
the view that purported weakness of spin-down has enabled the Sun to have sustained 
a high gravitational oblateness, although they do not tell us how. 

The situation is changed dramatically once it is recognized that the radiative en- 
ight be pervad ed b y a large- s cale m agnetic field. Mestel and Weiss (1987; see 
lll953l.l2012l'l and lCowlind lll976h ') have argued that a quite modest poloidal 
field (about 5 x 10~®T) is sufficient to connect the core to the convection zone in the 
lifetime of the Sun, and have advanced arguments to suggest that the actual field could 
be of order 10~^T. That field is substantially weaker than the only ab initio, albeit 
poor, estimate - of order SOT - I Goughlll990al ) t hat f know. The magnetic spin-down 
process was subsequently studied numerically bv ICharbonneau and MacGregoJ 1 19921 . 
ll99Sbl[T993al l with an idealized model having a rigidly rotating convection zone. That 
study confirmed that magnetic spin-down is efficient, but the model was too simplistic 
to explain why the radiative interior rotates uniformly, an issue to which I now turn. 


velo; 

also 


pe migt 

iMestell 


6 The steady laminar tachocline 


The angular velocity of the convection zone varies with latitude in an almost depth- 
independent manner, and interfaces with a nearly uniformly rotating radiative interior 
via the thin tachocline. How is that achieved ? The matter was first seriously considered 
by Spiegel and Zahn (1992, see also ISniegell ( 1972|H . They, like I here, did not discuss 
the cause of differential rotation of the convection zone: it is evidently the result of a 
balance between the angular-momentum-transporting Reynolds stress. Maxwell stress 
and advection by larg e-scale meridional flow, a matter which is reviewed briefly by 
iHanasoge et al.l 1 2015h in this volume. Instead, recognizing that the global equilibration 
timescale of convection is no doubt much shorter than the dynamical timescales of 
processes operating beneath (even if they are related to the solar cycle), Spiegel and 
Zahn considered the convection-zone shear to be given, and ignored any back-reaction 
of the tachocline dynamics on the convection zone. They then asked why the bulk of 
the fluid below rotates uniformly. 

Spiegel and Zahn first established that although the radiative interior is very highly 
stratified, thermal diffusion is sufficiently efficacious to permit baroclinically driven 
flow (induced in a manner similar to that in my spin-down discussion of the previous 
section) to advect the required angular momentum essentially throughout the Sun in its 
lifetime. Therefore, the (imposed) differential rotation of the convection zone has to be 
insulated from the radiative interior to allow the latter to rotate uniformly. To achieve 
that, Spiegel and Zahn invoked a tachocline pervaded by turbulence generated by the 
tachocline shear itself, suppressed vertically by the stable stratification and therefore 
being layerwise two-dimensional. Moreover, they implicitly assumed the turbulence 
to be horizontally isotropic, notwithstanding the fact that the rms vorticity of the 


p erhaps even in the face of the steep gradient of molecular weight dHuppert and Spiegell 
Il977fl which had previously been re garded as isolating the angular momentum of the energy¬ 
generating core llMestellll953l . Il957f l. 
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turbulence is likely to be comparable with, or perhaps even less than, the angular 
velocity of the Sun, a situation which one would naturally expect to lead to substantial 
anisotropy. Thus they introduced a (turbulent) viscous stress tending to force the 
rotational flow towards being horizontally shear-free, and were thus able to achieve 
a steady state with a thin tachocline abutting a uniformly rotating interior, having 
angular velocity 17c. In the tachocline itself there was a flow similar to that in the 
Ekman layer in Figure 5: towards the axis of rotation - i.e. poleward, on essentially 
horizontal surfaces - in an equatorial region where the tachocline rotation exceeds that 
of the interior, and equatorward in polar regions. Spiegel and Zahn estimated that 
the two regions meet at latitude 42°. That implied that 17c is essentially the value 
of the photospheric angular velocity at that latitude, namely 0.90l7e, where 17e is the 
photospheric angular velocity at the equator. 

Th e study stimulated a potentially interesting numerical simulation by iMieschl 
( 2003l l . who set out to determine, amongst other matters, whether the shear turbu¬ 
lence would actually be likely to induce the rigidity produced by a scalar viscosity. He 
generated turbulence by imposing a grid of point sources of (gravity waves) in a stably 
stratified fluid under gravity. The grid was forced to rotate rigidly, and therefore the 
resulting mean flow inevitably tended towards rigid rotation too, in just the same way 
that gravity waves generated by wind over mountains exerts a drag on the wind (in 
the frame of reference in which the mountains are stationary). The conclusion that the 
turbulence leads to rigid rotation was therefore unjustified. What would be much more 
revealing is to simulate a situation having the sources move with the flow, so that no 
external torque is applied to the fluid. 

Spiegel and Zahn’s study has no doubt introduced much of the pertinent dynam¬ 
ics of the tachocline. However, the details - indeed even the basic principle - have 
come under criticism. For example, it is unheard of that continuously maintained 
shear-generated turbulence ca n completely annul the shear that drives it. Indeed, an 
investigation by lEIlio t3 (Il997h demonstrated that even the magnitude of the turbu¬ 
lent stress generated would be insufficient for requirements. Moreover, in the observed 
natural environment, predominantly the Earth’s (rotating) atmosphere, layerwise two- 
dimensional turbulence leads to augmentation, rather than suppr ession, of larger-scal e 
shear, partly through angular-momentum transpor t by waves i Havnes et al. |[l99l|). 
Such considerations led iGough and Mclntvra (Il998l i to conclude that the angular ve¬ 
locity of the radiative interior could never by kept uniform by fluid-dynamical processes 
alone, and that the interior must necessarily be held rigid by a large-scale magnetic 
field, presumably primordial. They outlined a nonlinear dynamical balance in what 
they regarded as the simplest model of the tachocline; it again led to a meridional 
flow, both geometrically and dynamically similar to that inferred by Spiegel and Zahn, 
associated with which is downflow from the convection zone in both polar and equa¬ 
torial regions, and upflow between, near the latitude of zero shear, observed seismo- 
logically to be about 30° (see Figure 4), and implying that f7c — 0.93f7e (cf. Figure 
2). Advection of magnetic flux by the downflow counters the tendency for the field to 
expand by diffusion, yielding a steady-state balance which determines the thickness of 
the tachocline. Near the shear-free latitude the outflow might lift the primordial field 
into the convection zone, anchoring the angular velocity in the convection zone to that 
in the radiative interior, and p ossibly fuelling the m agnetohydrodynamical processes 
responsible for the solar cycle ( Bvington et al.ll2014l ~l. An obvious consequence of this 
picture is that in this region the shear might actually be quenched completely by the 
penetrating magnetic field, a feature that in principle is testable seismologically. It is 
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interesting to note that the latitude of this region is the same as that at which sunspots 
first appear at the start of a solar cycle. It would be an unlikely coincidence if that had 
no dynamical significance. 

Of course, the potential role of a large-scale magnetic field holding the radia- 
tive interior rigid had b een discussed before. I have already mentioned the w ork of 
iMestel and Wei^ ( 1987ll and icharbonneau and ^ MacGr egoi] (|l993al jjf93bl [l993ah in con¬ 
nection with spin-down. Additionally, Rudiger and Kitchatin^ (Il997l l imagined the 
presence of such a field, without careful regard to the direct action on it of the dif¬ 
ferentially rotation c onvection zone, and proposed the tachocline t o be just a simple 
Hartmann layer (e.g. lHartmamilll937l : ISherclifjl 19651 : lRobertslll967l l as an explanation 
of why it is so thin. What Gough and McIntyre argued is that the magnetic field with 
a poloidal component must necessarily be present in the radiative zone, and they pro¬ 
vided an albeit approximate analysis of how it is prevented from crossing (most of) the 
tachocline by gyroscopic pumping from the convection zone, thereby disenabling it from 
imparting on the radiative interior the latitudinal rotational shear in the convection 
zone. Some aspects of the analysis were subsequently supp orted by numerical simu¬ 
lation s carried out by iKitchatinov and Riidigen ( 2006h and iR.iidiger and Kitchatino-^ 
(l2007lb in particular the advection domination by the meridional flow, although the 
models differ in their details. Rudiger and Kitchatinov assumed a meridional flow 
much more rapid than the value derived by Gough and McIntyre, and produced a 
tachocline above an essentially uniformly rotating interior with a less intense field. 
An essential feature of both descriptions is the absence of Maxwell stress at the base 
of the convection zone. There are more complicat e d pictures in which t he tachocline 
is magnetized and turbulent (e.g. iMcIntvra l2007l : IPiamond et al.l 1200711 . a result of 
the nonlinear breakup of the flow resulting from the magne torotational instability 
( Velikh^ WMyChandjasekhm 1960|: Balbu s and Hawlevlll99lfl or the Tay ler instabil- 
i tv (iTayl CT 1957l~ 19731: Markev ar id Tavleif 1973h - see al^ lSpruitl (120021) : lArlt et all 
( 2007 '!: Kitchatinov and R.udigeiT i 2007 1: Riidiger and Kitchatinov ll2007ll - and these 
are perhaps more representative of reality. 

Baroclinic meridional flow in the tachocline is an inevitable consequence of the 
basic steady-state dynamics of any model in which Maxwell stresses play only a minor 
role within the body of the tachocline itself. Therefore a prediction of the latitude 
of the confluence of the poleward and equatorward flows provides an easily testable 
consequence of the theory, because it determines f2c, which has been measured seismo- 
logically. That latitude in the Spiegel-Zahn theory is obtained from a linear equation, 
and is well determined once the latitudinal dependence of the turbulent stress-strain 
relation is specified (Spiegel and Zahn assume it to be constant, without comment). 
By contrast, assuming a local linear stress-strain relation across the tachocline, with 
a coefficient of proportionality that is independent of latitude, yields Qc ~ 0.96f7e 
( Goughlll985l ~). The Gough-Mcintyre theory is fundamentally nonlinear, and has not 
yet been worked out in sufficient det ail for a prediction to be ma de. Some simpler toy 
models have been looked at too fe.g. iGaraud and Guervill^l2009l 'l. which at least shed 
some light on the potentially relevant processes: vertical shear, horizontal advection, 
geometry of Maxwell and Reynolds stresses. 

One of the issues that needs to be addressed in the Gough-McIntyr e picture is how 
the co nfigu ration could have been establishe d. Numerical simulations bv iBrun and ZahnI 
(I2OO6I I and Istrugarek et al.l 1 2011all2011bl l. for example, have failed to reproduce an 
appropriately confined interior magnetic field, leading the authors to conclude that the 
theory must be wrong. But the establishment depends critically on the dominant role 
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of field advection, which could not be reproduced by the simulations because, as in 
all such computations, realistically small diffusion coefficients could not be achieved. 
Greater success has been achieved in a programme led by Garaud, who was able to re¬ 
duce diffusion by restricting attention to axisymmetric configurations (a restriction also 
adopted originally in the discussion by Gough and McIntyre). One of the consequent 
difficulties that she found initially is that any field pe netrating the t achocline at the 
poles could not readily be advected equatorwards (e.g. lGaraudll20o3 ). leaving a large 
stress on the axis which locked the rotation of the convection zone and the radiative in- 
terior, even though an unlo cked steady-state solution was sub sequently shown to exist 
( [Wood and Mclntvrell2nilh . However, the most recent work I Acevedo-Arreguin et al.l 
I2OI3II has succeeded for the first time to reproduce a steady two-dimensional numer¬ 
ical configuration consistent with the Gough-Mcintyre theory, and provide a simple 
explanation of why previous numerical investigations had failed. 


The assumption of axisymmetry, adopted originally for superficial simplicity, is 
perhaps also an extremely inhibiting impediment to the success of many of the the 
simulations. Provided other completing influences cannot dominate, an axisymmetric 
magnetic field, for example, whose axis of symmetry is initially nearly but not exactly 
aligned with the rotation axis, could perhaps be advected away from the rotation axis 
by the equatorward meridional flow in the tachocline, because it is on the magnetic 
axis that the field strength is greatest. The field would be left, maybe, with its axis of 
symmetry intersecting the tachocline at the latitude of the poleward-equatorward con¬ 
fluence: the latitude of zero shear. There it penetrates into the convection zone, where 
perhaps it suppresses the tachocline shear in a non-zero range of latitude. A seismolog- 
ical test for such an outcome is currently being developed. Indeed, it is not wholly out 
of the question that a similar process is at least partially responsible f or the obliquit y 
of the orientation of large-scale magnetic fields in earlier-type stars ( Goughll2012bl l. 
More complicated asymmetric configurations can be envisaged. The outcome would 
be essentially a steady magnetic field configuration rotating with angular velocity f?c, 
which overall leads to a nonaxisymmetric Sun, even if the field itself, in its appropriate 
frame of reference, were axisymmetric. Could this be the root of an explanation for the 
existence of active longitudes? The asymmetry would rotate with the Sun, maintaining 
its phase. However, long-term ph ase stability of the a ctive longitudes is not obviously 
evident in the observations (e.g. iGvenge et al]|2014h . Nevertheless, other indicators 
could be more telling: it is interesting that the sector structure of the solar wind, 
which might mirror the magnetic asymmetry within the Sun, has maintained its phas e 
for at least four sunspot cycles prior to the mid seventies ( Svalgaard and Wilcox|[l975l f. 
Further, similar, analysis of solar-wind data up to the present day would therefore be 
very welcome, and could add (or subtract) credence to the hypothesis. 


Descriptions of the genre that I have discussed have all led to a ventilation timescale 
by the baroclinic meridional flow of order 10®y. That is long enough for thermal dif¬ 
fusion to establish a radiatively balanced thermal stratification, yet too short to be 
countered by microscopic diffusion, or gravitational settling of chemical species. There¬ 
fore the tachocline is chemically homogeneous with the convection zone, and thus has 
lower mean molecular mass than it would have had had gravitational settling been 
unopposed. Therefore the sound speed is higher than it would otherwise have been, 
creating the positive anomaly prominent in Figure 6. lElliott and GoughI 1 1999ll used 
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Fig. 6 Optimally localized av erages of relative differences betwee n the squared sound speed 
in the Sun and in Model S of IChristensen-Dalsgaard et al.l (Il996lh computed by M. Takata 
from MDI 360-day data and plotted against the centres x = f/R of the averaging kernels 
A{x; x) — which here resemble Gaussian functions - and are defined hy x = J xA^dx / j A^dx. 
The length of each horizontal bar is twice the spread s of the corresponding averaging kernel, 
defined as s = 12 J{x — x)^A^dx\ an averaging kernel A that is well represented by a Gaussian 
function of variance has spread approximately 1.7A ~ d.72FWHM\ were it to be a top- 
hat function, its spread would be the full width, which is why s is so defined. The vertical 
bars extend to ± one standard deviation of the errors, computed from the frequency errors 
quoted by the observers assuming them to be statistically independent with zero mean. The 
sharp positive anomaly centred at x ciii 0.67 immediately beneath the convection zone is no 
doubt the consequence of chemical homogenization of the tachocline with the convection zone. 
The outward decline in 5c^ in the convection zone is a result of having underestimated the 
seismic radius of the Sun. The convex variation abo ut x 0.15 provides a hint of there having 
been some large-scale meridional flow in the core (cf. lGough and Kosovichevlll99Qh . which may 
also be responsible for the low equatorially biassed angular velocity evident in Figure 2. The 
broad positive discrepancy in the radiative envelope is not understood. 


that property to calibrate the thickness A of the tachoclinJ^. yielding A ~ 2.0x 10 ^R, 
which is more precise (and smaller) than seismological measurements of shear, because 
hydrostatic stratification can be measured more precisely than rotation; however, the 
outcome depends on the reliability of the value of the diffusion coefficients required 
for evaluating the extent of gravitational settling, so the estimate may not be as ac¬ 
curate. Indeed, a yet unpublished investigation by Christensen-Dalsgaard and myself 
has established that the precise form of the tachocline anomaly cannot easily be repro¬ 
duced by standard spherically symmetrical solar-structure theory, suggesting, perhaps, 
the presen ce of a degree of aspher icity of the tachocline, in contradiction to an earlier 
finding of iBasu and Antial (l200lll . From an assessment of the horizontal balance of 
forces, it is inconceivable that the base of the tachocline is both steady and aspheri- 


Here I adopt the original iSpiegel and Zahi] (119921) definition of the tachocline: the gyro- 
scopically pumped shear layer confined to only the stable region beneath the convection zone, 
despite the etymology of the appellation. 
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cal, although the transition between the stably stratified tachocline and the unstably 
stratified convection zone, aside from the inevitable convective buffeting, could be. 

Finally, a word about the overall shear near the base of the convection zone. Most 
seismological investigations of the tachocline have - quite naturally, given its appella¬ 
tion - considered only the shear itself, and have tried to characterize its thickness by 
fitting to t he ill resolved inv ersions for fl a chosen functional form with a quantifiable 
width (e.g. lKosovichevlll996l l. Not surprisin gly, the values found tended to exceed that 
determined from the sound-speed anoma ly jKosovichevl IlOflbl : iBas ul ll997l : lAntia et al.l 
1 19981 : ICorbard et al.l[l998l : ICharbonneau et al.lll999l l. because. in accord with the orig¬ 
inal de finition, the true tach ocline shear exists only in the stably stratihed interior, al¬ 
though [COTbarJ^eT^^ ( 1999h suggested subsequently that A might be as low as 0.017?. 
Almost all theoretical studies to date have ignored the reaction of the convection zone 
to the tachocline shear, which must penetrate to some degree into the convection zone, 
especially at high latitudes where vortex stretching is the greatest. Of course, some 
upward penetration of the vertical shear must occur, and is indeed clearly visible in 
the seismological inferences (Figure 2), especially in the polar regions where vertical 
shear is resisted the most strongly by vortex s tretching. It i s that penetration that has 
led to the conclusion by K qsoyich^ 1 1996l l: lAntia et al.l (Il998h : ICharbonneau et al.l 


(I 1999 I I and iBasu and Antia (1200,811 that the tachocline is thicker at the poles. A recent 
detailed heli oseismic study o f the s hear, including its temporal variations, has been 
presented by lAntia and Basul 1 201lll . 


7 The temporally varying tachocline 

The description of the tachocline outlined in the previous section is dominated by 
essentially laminar dynamics, superposed upon which there might be some weak small- 
scale turbulence providing additional, diffusive, transport. Moreover, the bulk of the 
tachocline is usually considered to be free of magnetic field, except in the upwelling 
region of near-zero shear. But there are more complicated situations that have been 
considered, and which suffer temporal variation on a timescale much shorter than the 
10®-year ventilation time mentioned above. 

The most obvious time variation to consider is the buffeting by large-scale plumes 
from the convection zone. This is likely to cause the stability boundary - i.e. the sur¬ 
face across which the (local) convective stability changes - to undulate, in a manner 
similar to the tops (and bottoms) of terrestrial clouds (which, it should be pointed out, 
are not demarcated precisely by the visible vapour interface). Beneath that interface 
are tight ripples - gravity waves predominantly with timescales and horizontal length- 
scales comparable with those of the buffeting - because the stable stratification is very 
much more intense than the unstable stratihcation above, with vertical lengthscales 
very much less than the horizontal scales. The group velocity is very nearly horizontal, 
so the ripples penetrate only a very short vertical distance before they dissipate (e.g. 
Gough, 1977). The undulations are usually called overshooting, but they do not nec¬ 
essarily induce as much mixing as is often presumed for causing the region of nearly 
adiabatic stratification to penetrate downw ards and terminate in a yet sharper inter¬ 
face (e.g. K^ristensen-Dalsgaard et al.lboilh '). Instead they induce a smoothing of the 
horizontally averaged stratification (to which global seismology is sensitive). There is, 
however, a small degree of material mixing resulting from minor disruptions to the 
interface. The gravity waves themselves also contribute to material transport by (non- 
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linear) Taylor (fl953l. 1953) dispersion, thouRh on a mu ch smaller vertical scale (e.g. 
IPress and Rvbickil 198ll : iKnobloch and Merrvfieldlll992l ). It is likely that the pressure 
perturbations associated with the convective plumes are transmitted to the tachopause 
- the base of the tachocline - causing undulation in the boundary beneath. A conse¬ 
quence of the smoothing of the horizontally averaged stratification near the base of the 
convection zone is a reduction in the magnitude of the (negative) leading constant i n 
the asymptotic formula for the periods of high-order g modes of low degree ( Ellisl[l983 ) . 
which could be detectable should the modes be observed. At present there is no hope 
of a direct observation. However, indirect methods, such as seeking g-mode-induced 
undulations in p modes, the dynamics of which I refrain from discussing here, may one 
day be successful. 

The simple Gough-Mcintyre description of the tachocline is magnetic-field free, 
except where it interfaces with the primordial field in the interior, and also near the 
almost shear-free region where the tachocline flow is upwelling. There is undoubtedly 
also a magnetic field in the convection zone above, vacillating with the solar cycle with 
a characteristic 22-year period. If the temporally aver aged field were strictly zero, the 
field would not penetrate far into the tachocline fe.s. lGaraudlll993 ). which is why in 
their most elementary description the possible existence of such a field was ignored 
by Gough and McIntyre. But isn’t it more likely that there is a randomness in the 
vacillation that leaves a nonzero residue? The residual field crossing the tachocline 
would be wound up by the tachoclin e shear. The situation would be ripe for Tay ler- 
type instability, as lSprui3 ( 19991 . [20021) has advocated (see also [Diamond et al.ll2007l ). It 
would lead to additional transport by the ensuing turbulence, and it would modify the 
relation between the tachocline thickness and the large-scale interior magnetic field 
( Mclntvr^ l2007l ). I think it is fair to say that this very complex subject is still ill 
understood. It remains an area of active research, attracting not only those who regard 
themselves as solar physicists. Fo r more inform a tion I refer the reader especially to the 
book on the tachocline edited bv lHuehes et al.l ( 2n07h . and to subsequent publications 
by the contributing authors. 

It has even been proposed that the residual vacillating field leaking from the con¬ 
vection zone is alone sufficient to render the rotation of the radiative interior uniform 
(Forgacs-Dajke and Petrovay, 2000, 2002), without recourse to requiring the rigidity of 
a large-scale field. The idea seems to demand that the leaking field that is generated 
by dynamic action in the differentially rotating convection zone somehow takes on a 
rigid configuration. What must surely be the case is that the field takes on at least 
part of the convection-zone shear, which is then transmitted through the tachocline. 
So far as I am aware, not even a highly idealized model of the pertinent dynamics has 
been fully investigated. However, if one imagines a toy model in which a periodically 
oscillating source of field with zero mean generated in a convection zone with high, 
turbulent, scalar magnetic diffusivity rotating at the observed differential rate, and 
diffusing into the radiative interior below, where the diffusivity is very much lower, 
then the temporally averaged Maxwell stress on the interior does not vanish, and is 
such as to generate differential rotation in the radiative envelope of the same form, 
although not necessarily with the same amplitude, as that in t he convection zone . 
Solar-cycle-related time dependence is also not out of the question. lAntia et al.l ( 20131 ) 
have embarked upon a programme to study potential dynamical consequences of the 
seismologically determined variations of the angular velocity in the convection zone. 
It is the only truly dynamical study that stems directly from helioseismological infer¬ 
ence. The magnetic field required to produce the seismically observed angular-velocity 
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variations was obtained under the assumption that the only azimuthal force is solely 
a Maxwell stress, and in this initial study only the azimuthal balance of forces was 
considered. The analysis is therefore far from complete. However, a magnetic field was 
obtained which oscillates not only in the body of the convection zone but also in the 
tachocline. That result therefore reinforces the opinion that the early, steady, tachocline 
studies were grossly oversimplified. 

I conclude this section by mentioning the so-called tachocline oscillation, an appar¬ 
ently wave-like oscillation in the angular velocity near the equatorial plane immediately 
above and below the base of t he convection zone , with a period of about 1.3 years. It was 
discovered seismologically bv iHowe et al ] (l2000l'l . with opposite phase in the convection 
zone from that at the top of the radiative interior, thus exhibiting a maximum of the 
shear amplitude in the tachocline. Soon afterwards the oscillation almost disappeared 
( Howe et al l lsoTih; the ey e of faith of an imaginative observer might perceive hints of 
its return isee lHo\rol2009ll . although the evidence is not statistically significant in any 
plausible sense. What is the restoring force? Howe et al. suggested a magnetic field, 
which is quite plausible. Indeed, there is (unpublished) evidence that the oscillation 
extends more deeply than has been reported in the literature, with a further maximum 
in the shear at about r — 0.55R, but with rather lower amplitude, as one would expect 
if the field does not increase with depth faster than y/]}. It is interesting, although 
perhaps merely coincidental, that the intensity of the vertical component of the field 
required to produce such an oscillation would need to be about IT, which is similar 
to, yet rather less than, the rough ab initio estimate of the global field in the radiative 
zone 0 Howe and her colleagues (personal communication) consider the evidence for 
the deep penetration of the oscillation to be too insecure for them to report, but in 
trying to understand the workings of the Sun it is useful for theorists at least to bear in 
mind the possibility. What might drive the oscillation? One possibility is its influence 
on the anisotropic dissipation of gravity waves, a process which drives the terrestrial 
QBO, to which I turn my attention later. 


8 Excitation of seismic modes 


It is generally agreed that in the Sun acoustic oscillation modes are intrinsically sta¬ 
ble, their energy being absorbed into the background configuration of the star by the 
combined action of appropriately phased heat and momentum transport by convection 
and, to a lesser degree, radiative heat transfer in the convectively stably stratified in¬ 
terior. The modes are generated stochastically by the turbulence in the upper layers 
of the convection zone, by both random impulses from the turbulent motion and by 
the associated fluctuations in buoyancy. There has been a sequence of more-and-more 
sophisticated analyses of the processes involved (e.g. Stein|[l96^ Goldreicl^_and_Keel^ 


19(^: fBalmforth and Goughlll990l : lBalmforthlll99a: Bej^re^eTa^^OO^ Cha 2 linet_^ 


2005). They derive initially from the work of Batchelor (I 1953 II and iLighthilll (I 1952 I . 
1954l ~l. who considered respectively the stochastic excitation of a (general) simple har¬ 
monic oscillator and the mechan i sm of turbul e nt acoustic wave generation by mo- 
mentum transfer - also of iGoughl ( 19651 . Il977al l: lUnnol (Tio^); IXiond 1 1989ll : iGabriell 
(Il99fih . who addressed the role of convection in determining the linear growth and 


and similar also to field in t ensiti es commonly discussed elsewhere in connection 
with the tachocline (e.g . [SBrui ^ |2002|: [Gilman a nd Callvl 1200711 . although others (e.g. 
iGough and McIntvrdfTooS:|Kitchatino^an^^nidige^l2007r) have entertained weaker fields. 
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damping rates of stellar oscillations (see also lPunret et aDboOiil . boOfj . l2009li . The ab¬ 
solute intensity of the nonlinear excitation is extremely difficult to quantify, because 
it is extremely sensitive to details of the turbulence, which are not well defined by 
the mixing-length fo rmulations that are adopted to represent the solar convection zone 
I Goughl[l977al . l2002h . Therefore the absolute amplitudes of the seismic waves are very 
uncertain. Nevertheless, an appropriate combination of the uncertain factors in the 
theory can be adjusted at least to harmonize with seismic observations; that leads to 
a functional form of amplitude against frequency - a true prediction of the theory 
- wh ich is in reasonable accord with observation I Houdek et al.l [l99fll : IChanlin et al.l 
l200fil l. It is interesting, and perhaps somewhat surprising, that, once the adjustment 
has been made, the remaining weakest part of the theory, at least when applied to the 
Sun, then appears to be the prediction of the linear damping rates of the modes; that 
inference was obtained by replacing the theoretical damping rates by o bserved acous¬ 
tic p ower-spectrum line widths, and finding substantial amelioration ( Chaplin et al.l 

I2OO5II . 

Further work on the convection-pulsation interaction is therefore called for, not 
merely to make minor improvements to derivations of mode amplitudes in the Sun, 
but, more importantly, to improve our understanding of the underlying dynamics, both 
for its own sake and for application to other stars for which the gross solar adjustment 
of the theory may no longer be appropriate. 

Observations of solar-like observations in other stars add new data with which to 
compare the theory. Particularly useful are stars rather different from the Sun. An 
observation of the mode amplitud es in the (coolish) more e volved star ^ Hydrae, for 
example, has been well predicted ( Houdek_and_^ugh| 2003 1 . although subsequent di¬ 


rect estimates of the damping rates (IStello et al.ll2006li ' were not in accord with the 
theory, suggesting, perhaps, the presence of cancelling errors in the theoretical growth 
and excitation rates. Before jumping to that conclusion, however, one must recognize 
that the damping rates were obtained from the observations by equating them with 
the line widths in the power spectrum, which is a correct procedure only if it is known 
that the background sta te of the star is truly invariant over the duration of the obser¬ 
vations I Batchel^ll953ll : otherwise non-dissipative phase (frequency) wandering can 
broaden the spectral lines and thereby falsify the conclusion. Indeed, subsequent analy¬ 
sis of high-quality d ata from the space missions CoRoT and Kepler ( Huber et al ] l20inl : 
iBaudin et al.ll201ll 'l yielded values in accord with the theoretical predictions. 

The situation is much w orse with the hot sta r Procyon, whose oscillation am plitudes 
are grossly overestimated I Houdek et al ] ll99fll : [Matthews et al.ll2004bll2004al l. In this 
case there must be something fundamentally wrong with how one describes either how 
the turbulence properties scale with variations of stellar structure, or, worse still, the 
manner in which they combine to quantify the excitation; alternatively, there might be 
some missing ingredient in the calculation of the damping rates, such as wave scattering 
or some process unaccounted for in the time-dependent convection theory. 


9 On the dynamics of gravity modes 

In an attempt to resolve the solar neutrino problem before helioseismologv. IPilke and GoughI 
(Il972l) suggested that intermittent nonlinear breakdown of a grave unstable grav¬ 
ity mode would trigger finite-amplitude convection that almost completely mixes the 
core. That suggestion is now known not to be correct, because the form of the de- 
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pression of the sound-speed in the central regions of the Sun, evident in Figure 1, 
incont rovertibly demonstrates the e xistence of a substantial gradient of helium abun¬ 
dance ( Gough and Kosovichevl[l990[l . Moreover, an important analysis bv IPziembowskil 
( 1 9821 1 of triad coupling to a resonating pair of otherwise damped daughter modes 
suggested that it is unlikely that the principal unstable mode would ever achieve an 
amplitude high enough to overcome the barrier to the nonlinear onset of convection. 
Nevertheless, the dynamical processes that were thought to operate are not necessarily 
irrelevant today, and can perhaps be profitably rehearsed. The driving of grave solar g 
modes is prin cipally by the mo dulation of the nuclear reactions, a mechanism under¬ 
stood first bv lEddingtonI ( 192fJ) . and now called the e mechanism. Extreme sensitivity 
of the nuclear reaction rates to temperature causes (local) temperature excesses to be 
augmented in the face of (global) losses by radiative diffusion. Key to the process in cool 
stars like the Sun, which are powered by the relatively insensitive p-p reaction chain, 
is the destruction of nuclear balance by the oscillations, which vary on the timescale of 
an hour, much shorter than the 10^-year nuclear equilibration time, exposing the varia¬ 
tion of the energy ge nerated by the temperature-sensitive ^He-^He, and the ^He-^He, 
nuclear reactions fcf. ILedoux and Sauvenier-GofMll95fll l. Also important is that the 
g-mode-depressing convection zone is deep compared with those in hotter stars, pre¬ 
venting the mode from achieving so great an amplitude in the outer stellar envelope as 
to suffer sufficient damping to overc o me the nuclear dr i ving. C onsistent calculatio ns by 
Christensen-Dalsgaard et al ] lll974h . [shibahashT et al ] (Il97,’il 'l , iBourv et al.l (ll97!Th and 
Said ( 198ol l confirmed the original estimate that some g modes were indeed unstable. 


The analysis bv IPziembowskil ( 1982l l is particularly interesting dynamically, because 
the outcome is superficially counterintuitive. It is well known that a pair of daughter 
gravity modes can sap energy from their parent, provided that they resonate precisely 
enough. What is less well appreciated are the conditions required for that process 
to be severely debilitating. When I was a student at the Program of Geophysical 
Fluid Dynamics (GFD) at the Woods Hole Oceanographic Institution in 1965 my 
attention was caught by the daily water being heated for coffee in an aluminium pan 
on a fiat electric hob. The bottom of the pan was rounded with age, and the pan 
rocked on the heated plate. If the quantity of water was right, there was a resonance 
between the rocking frequency and the frequency of surface gravity modes in the water, 
and the differently located bursts of heat through the bottom of the pan as different 
parts of it touched the hot plate caused the oscillation to be sustained. Eddington’s e 
mechanism was operating at home. The oscillation amplitude was not tiny, yet surely 
the intrinsic dissipation - measured by the inverse Reynolds number, for example - 
was much greater than it is in the Sun. So shouldn’t gravity modes in the Sun attain 
even greater amplitudes? The answer is negative: oddly enough, it is actually because 
the propensity for dissipation in the GFD pot was so high that the oscillation could 
be sustained. To understand that one must first realize that to be an effective energy 
drain the difference between the frequencies of the two daughters must match the 
frequency of the parent within an inverse intrinsic growth time, a condition which in 
general becomes harder and harder as the thermal diffusion coefficient, and hence the 
nonadiabaticity, diminishes. This requires finding daughter pairs amongst modes with 
higher and higher characteristic wave numbers, which actually leads to higher and 
higher damping rates, notwithstanding the smaller diffusion coefficient, therefore lower 
and lower filial amplitudes and consequently an even lesser and lesser ability of the 
daughters to extract energy from their parent. Only if the dissipation rate (diffusion 
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coefficient) is low enough can the daughters attain high-enough amplitudes for them 
to be able to limit the growth of the parent effectively. 

Dziembowski assumed the background stratification of the Sun would not be chang¬ 
ing with time. He found that it was not possible always to find daughter pairs suffi¬ 
ciently close to resonance, yet with adequately low damping rates to permit them to 
grow to amplitudes high enough to severely inhibit the parent. So he carried out a 
statistical calculation to estimate the expected amplitude of the parent, and found the 
amplitude to be a mere 10 cm s~^ in the early stages of the Sun’s evolution on the 
main sequence I Dziembowskil [TfiS.ll 'l . That is surely insufficient to cause serious dis¬ 
ruption to the structure of the core. Nevertheless, it appeared to permit accidentally 
much higher amplitudes. What is the possible outcome? That question is answerable 
once it is realized that the Sun is evolving, and that resonance with particular daugh¬ 
ters cannot be maintained. In fact the passage from one resonating pair to another as 
the Sun evolves typically occurs on a timescale shorter than the growth time of the 
daughters. Moreover, it turns out that strict formal phase coherence over the time it 
takes for the parent to reach its limiting amplitude is not required: the parent’s phase 
drifts as the Sun evolves, and as new energy-sapping daughters are encountered their 
phases are established by th e triad interaction itself. Therefore, as reported briefly by 
Ijordinson and Goughl ( 2000l l. the limiting amplitude of the parent, as its frequency 
drifts amongst the family of resonances, is simply the expectation value calculated by 
Dziembowski. A simplifying assumption of the original calculation was that interac¬ 
tion with only a single daughter pair was considered; in practice, however, there are 
interactions at any time with an entire hierarchy, limiting the amplitude of the parent 
yet more severely. However, in general it is the gravest pair who sap by far the most 
energy, and it is likely that the extra loss to all the other daughter pairs is actually 
quite small, leaving Dziembowski’s estimate more-or-less intact. 

It appears from this discussion that the structure of the Sun, even the temperature- 
sensitive ®B neutrino production, is unlikely to have been modified seriously by grav¬ 
ity modes. Indeed, an end to this matter has been called by many astrophysicists 
who have grappled with neutrinos. However, anyone interested in dynamics might still 
ask whether there is any small remnant conseq uence. There has for a long time (cf. 
iGough and Kosovichevl[l990l : iGough et aklligOa l been a hint, albeit rather weak, that 
some degree of material redistribution, possibly laminar, has taken place in the in¬ 
homogeneous energy-generating core: this is evinced by the generally negative slope 
of !(? inwards of r ~ 0.2/?, visible in Figure 6, which appears to imply a radial 
gradient in the spherically averaged Y that is somewhat lower than theoretical expec- 
tation, and which , both in the very early days of helioseismology and also more recently 
( Basu et al.ll2001^ b has been recognized as being suggestive of mild mixing. Moreover, 
a general redistribution of angular momentum within the core from an initially uni¬ 
formly rotating state could make the core appear to be rotating more slowly, as is 
suggested by Figure 2, because on the whole the angular velocity near the edge of the 
core, especially near the equatorial plane, where rotational splitting is most prone to 
be affected, is likely to be diminished. Could these hints be manifestations of a soliton 
in the core, controlled either by gravity-wave-like dynamics, or by some more slowly 
operating d ynamics that m ight be influencing reported variations in the neutrino lumi¬ 
nosity (e.g. Sturrocl3l2008l l ? A possibly simpler, but less interesting, solution might be 
that mixing in the early stages of the main-sequence evolution was much more exten¬ 
sive than it is now, because the angular velocity (and especially its deceleration) was 
much greater then than it is today, possibly causing adequate baroclinicity to drive a 
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laminar circulation through the core when the impedance by chemical inhomogeneity 
was quite small. 

Could there have been significant Taylor dispersion by daughter g modes? These are 
matters that appear not to have been addressed. One might note that the grave g modes 
have similar amplitudes both in the core and near the surface, and that Dziembowski’s 
10 cm s~^ estimate is high enough to be observed at the surface today. However, the 
structure of the Sun today differs substantially from its early m ain-sequence state , 
and some calculations suggest that no g mode is now unstable I Bourv et al.l HOT!^ : 
IShibahashi et allllOTlih . although the matter is a delicate balance between driving 
and damping process es, leaving open the possibili ty that some g modes are today 
intrinsically unstable dNoels et al.lll97a : ISai 3[l98fl). There have been various claims 
of direct detection, but these are not universally accepted. I think it is fair to say 
that there is yet no convi ncing evidence for their existence at observable amplitudes 
( Appourchaux et al.ll201(]r) . Finding gravity mo des would be extremely importan t for 
seismology; as has long been appreciated fe.g. lOoughlllQTsI : ICabriel et al.fll99!Tll . the 
additional seismological resolving power that the g modes would provide would be quite 
substantial. Are there modes that are unstable, or are they all intrinsically stable and 
excited stochastically as pressure fluctuations associated with the undulating tachocline 
or by the same pressure fluctuations that excite acoustic mo des in the outer layers of th e 
convection zone? Estimates of the latter are summarized bv lAnpourchaux et al ] (120101) . 
They vary quite widely, although amplitudes are typically lower than the observational 
claims. A great deal of uncertainty remains. 


10 Gravity-wave transport 


In an early d iscussion of gravity-wave transport between core and convection zone 
( Gough|[l977al ) it was estimated that grave modes of the kind that might be self-excited 
in the core to amplitudes comparable with those discussed in the previous section 
are too weak to carry a substantial amount of heat, and too weak also to transfer a 
substantial proportion of the Sun’s angular momentum in its lifetime. Moreover, stable 
gravity waves that resonate in space and time with the undulating lower boundary 
of the convection zone pene trate only ab out 1 km into the radiative interior. A more 
recent estimate of the latter I Goughl[2002l ) suggests that penetration is about 100 times 
deeper, but that is still too shallow to have a material impact on the overall structure. 

It was estimated that modes generated off-resonance would have too low an amplitude 
to transport a substantial amount of angular momentum. It was also pointed out that 
any group of gravity waves that do penetrate must dissipate asymmetrically, and are 
expected to enha nce, rather t han suppress, differential rotation. 

Subsequently [p^ 3 (Il98ll) and lCarcia Lonez and Soruitl (Il99lh estimated the wave 
amplitudes by balancing inertially generated pressure fluctuations near the base of 
the convection zone against corresponding resonant gravity-wave fluctuations in the 
radiat ive zone beneath. As with the estimates of acoustic wave generation (e.g. lOoughl 
Il977al ). there is substantial uncertainty arising from uncertainties in the properties of 
the convective turbulence. 

An important consequence of the enhancement of rotational she ar by gravity-wave 
transp ort was illustrated with a pioneering experiment carried out bv IPIumb and McEwarJ 
( 19781 ). who demonstrated that the necessarily anisotropic wave dissipation can induce 
global oscillations of the shear, with characteristic period inversely proportional to the 
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fourth power of the amplitude of the gravity waves (and proportional to the sixth power 
of the horizontal phase spe ed). The experiment was c arried out to confirm a theory of 
Ihindzen and HoltonI ( 1968h and iHolton and LindzenI ( 1972h explaining the mechanism 
of the terrestrial QBO. It is of some interest to summarize the basic process here, 
because there is misinformation in the literature. 

Low-amplitude gravity waves generated with a broad spectrum of frequencies at an 
essentially horizontal surface propagate upwards (in the case of the Earth’s atmosphere) 
into the sheared stratosphere. The large-scale flow in the stratosphere is essentially hor¬ 
izontal, its variation with time and horizontal distance being much weaker than the 
frequency and horizontal wavenumber of the waves (and was therefore ignored in dis¬ 
cussing the linearized local wave dynamics). Frequency and horizontal wavenumber in 
an inertial frame of reference are therefore conserved. Given that the gravity-wave dis¬ 
persion relation is approximately = N‘^k^/{ky -l-fc^), where a; is the wave frequency 
in the local frame of the fluid, N is the buoyancy (Brunt-Vai'sala) frequency of the 
stratihcation, and and fev are respectively the horizontal and vertical components 
of the wave number, it follows that the horizontal component of the group velocity, 
and the (angular) momentum flux, are in the same direction as the horizontal phase 
velocity Also prograde waves, having negatively Doppler-shifted frequencies cu, have 
higher fev than do the corresponding retrograde waves, and therefore dissipate faster, 
causing the momentum that is imparted on the mean flow to be preferentially in the 
prograde direction. If prograde and retrograde waves are excited to the same ampli¬ 
tude, the shear is therefore at hrst enhanced, but further from the excitation layer, 
where the remnant prograde waves have lower amplitude as a result of earlier higher 
dissipation, the direction of the wave stress acting on the mean flow is reversed. 

An initially static atmosphere into which gravity waves are propagating is (then) 
susceptible to shear enhancement. If initially there is no mean shear, a random fluctu¬ 
ation in the symmetry of the wave excitation will first produce an initial weak shear 
which is then amplified by the process described above. The increased shear near the 
excitation layer eventually loses momentum to the excitation layer, and the whole pat¬ 
tern descends, producing a temporal oscillation in the mean horizontal flow at any 
given altitude. Evidently the amplitude and period of that oscillation must depend 
on the excitation rate and the dissipation rate, the latter, in view of the dispersion 
relation, depending on the horizontal phase velocity. The existence of the oscillation 
requires the fluid not to be too diffusive, for otherwise the oscillation is damped. 

Is there an analogous phenomenon in the Sun? Lik ely amplitude estimates lead to 
a characteristic period of about a month I Goughl[2002h . in a very thin layer within the 
tachocline into which the waves penetrate before being dissipated. It is most unlikely 
that any significance can be attributed to the coincidence of that period with the 
characteristic large-eddy turnover time at the base of the convection zone or with 
the rotation period of the Sun. In fact there may be no such coincidence, because 
the amplitudes of the waves generated, and hence the period of the oscillation in 
the shear generated, are so uncertain. A characteristi c period even a hundred times 
longer is not out of the question. Indeed. fOoughl 1 1998h wondered whether Plumb and 
McEwan’s (non-magnetic) mechanism could be responsible for the 22-year solar cycle, 
leaving the magnetic field to play only a secondary role in the dynamics. The process 
could also induce meridional flow l Eritts et aLlIlQQa i. which might conceivably provide 
a significant transport mechanism for chemical species, and possibly explain the Li 


The vertical component of the group velocity is directed oppositely to the phase velocity. 
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deficiency in the S un’s photosphere. An example of a rather different claim comes from 
an investigation bv iKumar et alJ 1 1999h . who estimated an enormously greater gravity- 
wave amplitude for waves that penetrate more deepl-J^. They too commented that 
Plumb and McEwan’s (non-magnetic) mechanism might be responsible for the 22-year 
solar cycle. In a subsequent investigation, Talon, Kumar and Zahn (2002) argued that 
shear instability would soften the oscillations. The effect of the instability was modelled 
with a scalar viscosity which increased the QBO-like oscillation period to about 300y, 
and the spectrum of the waves had enough power at low order to permit angular 
momentum to be redistributed throughout the radiative interior, again supposed to 
establish an almost uniform rotation profile. How the waves w ould suppress the imprint 
of the differential rotation in the convection zone discussed bv ISpiegel and ZahnI ( 19921 ') 
remained unexplained. 

The simplest discussions of this process are based on the action of an almost 
monochromatic group of waves (e.g. ikim and MacGregorl I2OO1I . l2003h . A more real- 
istic broad spectrum of waves such as would be expected in the Earth’s atmosphere 
{ Baldwin_et_ay[2001 ) and the Sun can lead to much more complicate d behaviour, 
Mathis_^^_aL|j 200^}_2ortend. From a ser ies of numerical simulations, Ro gers et al.l 


( 20081 1. [Rogers and MacGregor! ( 20ld . l201lf) and iMacGregor and Rogers ll201lh have 
concluded that simple linearized wave dynamics is inadequate for computing the mean 
Reynolds stress because (nonlinear) wave-wave interactions are important at the am¬ 
plitudes they expected to be generated near the lower boundary of the solar convection 
zone. They argued that oscillatory shear flow of the QBO type is muc h less like ly to 
be produce d - although rigid rotation is not to be expected either, as lMcIntvr'3 (e.g. 
1 19941 . 1201^ has often stressed. One might wonder whether in the numerical simulations 
excessive diffusion require s unrealistically large wave am plitudes, as is possibly the case 
also of convective motion i Toomre and Thompson 2015h. an d that the quasi-linear dy- 


Phimb and McEwa: 


namics of the kind described by 
await with interest the further, more realistic, simulations that are promised. 


3^ 1978ll really is a good guide. We 
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